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Using a sample of 12.2 million r-lepton pairs produced by e + e _ annihilation at yTs ~ 10.6 GeV and collected by the CLEO detector, 
we search for and set limits on CP violation in r-lepton decays. For each event, we require that both r-leptons decay via the mode 
r —¥ ttit !/. The search is performed within the context of a multi-Higgs Doublet Model and the imaginary part of the coupling 
constant parameterizing the non-Standard Model diagram leading to CP violation is constrained to be —0.046 < S(A) < 0.022 at 
90% CL. The novel search technique is of general utility. 



1 Introduction 

CP violation has long been observed in the kaon sys- 
tem and is expected soon to be seen unambiguously in 
the -B-meson system. Searches for CP violation in lep- 
tonic decays have been less extensive. Although CP non- 
conservation in such decays isJbrbidden in the Standard 
Model (SM), many extensions!!! to theJSM incorporate it. 
Multi-Higgs Doublet models (MHDM)Era are perhaps the 
best known. This report describes a novel search tech- 
nique for CP violation performed with tau decays and for 
a specific MHDM. The technique is general and no fea- 
ture of it restricts it to our particular choice of MHDM 
or final state. 

With the CLEO detector, we study e + e _ — > t + t~ 
events produced at CESR at a center-of-mass energy in 
the vicinity of the T(45) resonance. We select only 
those events where each r lepton decays in the mode 
t — ► TT7r ^ r . Interference between the SM process medi- 
ated by VF-boson exchange and a non-SM process medi- 
ated by scalar boson exchange could lead to CP violating 
effects. For each of our events, we construct a CP sen- 
sitive variable for which a non-zero value is evidence of 
CP violation. 

CP violation is essentially an interference effect be- 
tween two (or more) transition amplitudes that have a 
relative phase difference for both the strong and weak 
phases. The relative strong phase remains invariant 
under CP conjugation while the relative weak phase 
changes sign. Figure shows the basic idea, where S and 
</> are the relative strong and weak phases, respectively, 
and the overall transition amplitude A — A\ + ^e^e* 5 
for a process leads to its corresponding probability den- 
sity P: 
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cess, leading to difference in probability densities between 
a particular transition and its CP conjugate. 




A 2 e'^e'* 5 




Note that the last term contains the CP-odd factor 
sin 5 sin 0, which changes sign for the CP-conjugate pro- 



Figure 1: Generic Feynman diagrams for two transition amplitudes 
that contain CP-even and CP-odd relative phases. 



2 Model Choice 

To make the foregoing generalities specific and to illus- 
trate the overall search technique, we select a specific 
MHDM model: a 3 Higgs Doublet modelcl. The relevant 
Feynman diagrams for this model describing r — > irir is T 
are shown in figure @. The SM diagram mediated by a W 
boson is shown on the left while the non-SM diagram me- 
diated by a charged Higgs is shown on the right. For the 
SM diagram, the vector form factor that characterizes 
the final state strong interaction among quarks contains 
a CP-even phase. The non-SM diagram contains a scalar 
form factor f s , also with a CP-even strong phase. The 
choice of f s is not unique so we consider three choices: 
f s = l,f s = BW[a (980)], and f s = BW[a (1450)]. The 
term BW specifies a Breit-Wigner shape and ao indicates 
a particular intermediate meson resonance. 

In this model the second amplitude contains an over- 
all complex coupling A between the charged Higgs H and 
the final state, where A is a function of the H coupling 
to quarks and leptons. Because A is complex, it intro- 
duces a CP odd phase. Since we have two interfering 
diagrams with different relative strong and weak phases, 
CP violation is possible. 
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Figure 2: Feynman diagrams for a) the SM amplitude mediated 
by W exchange and b) the non-SM amplitude mediated by scalar 
exchange. 



3 Choice of Optimal CP-sensitive Variable 

In general, we can write the r production and decay prob- 
ability density P as the sum of two terms, one invariant to 
CP conjugation of the underlying Feynman amplitudes, 
Pevem and another that changes under CP conjugation: 
P = P even + Podd- Both P even and P od d are functions of 
the form factors, the complex coupling A and kinematic 
quantities. To search for CP-violation we select the CP- 
sensitive (i.e. CP odd) observable with the greatest sta- 
tistical significance and compute it for each event. Such 
a variablet3i3 has been constructed with different searches 
in mind by other workers: 
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£ is a complicated function, and although it is max- 
imally sensitive to CP violation, has no immediate intu- 
itive appeal. Averaged over the. data set, (£) ^ implies 
CP violation. It can be shown that 



(£) = Im(A) J |add. dLips 
= ciA + c 3 A 3 + c 5 A 5 + 



(3) 



where A = Im(A) and the a are constants. 

The Monte Carlo £ distribution is shown in figure pi 
for the two limits of the model. The left-hand side is for 
no CP violation and the right-hand side is for the 'maxi- 
mal' case where Jm(A) has a value equal to the SM cou- 
pling. The three distributions in each plot correspond to 
the three choices of the scalar form factor and the struc- 
ture in each distribution is due to the resonant structure 
of the vector and scalar form factors. A non-zero value 
of (£) for the CP violation case is apparent. 

4 Experimental Details 

4-1 Event Selection and Background Estimates 

To calculate a value of £ for each event it is necessary to 
know the direction of the final state neutrinos to deter- 
mine the direction of the final state r-leptons. For those 
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Figure 3: Monte Carlo distributions for the optimal CP-sensitive £ 
for the case a) of no CP violation and b) of 'maximal' CP violation. 
See text for discussion. 



events where each r decays semi-hadronically, energy and 
momentum conservation can be used to. determine each 
r direction up to a two- fold ambiguity a. Since we have 
no way of knowing the correct solution for an event, we 
sum the £ distributions corresponding to each solution 
and then average the result to search for an asymme- 
try. We have checked with Monte Carlo events that this 
process does not bias our result for the case of no CP 
violation. For the case of CP violation, we use a special 
Monte Carlo calibration procedure described below. 

We use a data set that corresponds to 13.3 fb _1 of to- 
tal integrated luminosity and contains 12.2 million t + t~ 
pairs produced from e + e~ collisions at a center-of-mass 
energy near or at the T(AS). We select events for which 
each r decays via the mode r — -» -kti v, selected for its 
large branching fraction (~ 25%)) and high selection ef- 
ficiency (10%). We verify with Monte Carlo events that 
our selection criteria do not introduce a bias in the £ 
distribution. 

The major source of background are r-pair events 
where one of the r's decays to a pv final state, the other 
r decays to a itit v final state, and one of the 7r°'s is not 
reconstructed. Using Monte Carlo events, we estimate 
this contamination level at 5.2% of the selected data. 
The overall background contamination from r decays is 
estimated at 10%. We estimate that events of the form 
e~e + — * qq (q = u,d.c, s quarks), e~e + — » T(4S) — > BB, 
and e~e + — ► e~e + 77 contribute a background contami- 
nation of less than 0.1%. 

4-2 Calibration 

From Eq. (g) we see that to lowest order (£) = ciA. To 
determine or limit A, we need to measure (£) and deter- 
mine the calibration constant c\. This is done in a two 
step process. First, we determine the range over which 
(£) is a linear function of A, using Monte Carlo events 
and assuming perfect detector response. Figure shows 
(£) v. A for the three choices of the scalar form factor. 



0.250 




3350301-006 


I i i i i I i i i i I i 


I | I I I I | 




A f s = 1 






_ □ fs = BW[aO(980)] 


/ 




• fs=BW[aO(1450)] j 


a 


0.125 




j*^*'*- 


A 

jj/> 

V 




- 


-0.125 




- 


-0.250 


I , , , , I , , , , I , , 


,!,,,,! 



-1.0 



-0.5 







0.5 



1.0 



Figure 4: (£} v. A for each of the three choices for the scalar form 
factor: a) f s = 1, b) f a = BW[a (980)], and c) f s = BW[a (980)]. 



After this linear range is determined, we choose 5 values 
of A within that range, and with full detector simulation 
compute (£) for each of the 5 A values. This produces a 
set of 5 (A, (£)) points to which a straight line is fit. The 
resulting slope is then c±. This entire process is done for 
each of the 3 scalar form factors. 

5 Results 

5.1 Observed £ Distributions 

Using the calibration constant, the £ distribution for data 
and Monte Carlo is shown in figure H for each of the scalar 
form factors. The data are squares and the line is SM 
Monte Carlo. The corresponding measurements for (£) 
and the 90% CL limits on A, the imaginary part of the 
non-SM coupling, are shown in table flL 

Table 1: Limits on Im(A) for various form factors / s . 



fs 


(0 x io- ;i 


Jm(A) @ 90% CL 


4 = 1 

is = BW[a (980)] 
f s = BW[a (1450)] 


-0.8 ±1.4 

-0.1 ±0.4 

0.1 ±1.2 


-0.046 < A < .022 
-0.008 < A < .006 
-0.019 < A < .021 



possible difference between track reconstruction efficien- 
cies for 7r~ and n + as a function of pion momentum are 
studied. The momentum distribution for charged pions 
in the reactions r — ► it^h^Vt is plotted in figure |6|a. 
The corresponding ratio of these distributions is shown 
in figure |6jb and is seen to be consistent with 1. Varying 
the slope of this ratio by ±lcr leaves c\ unchanged but 
does change A by ±0.003. We take this as a measure 
of the systematic error due to differences in the charged 
pion track reconstruction efficiencies. 

Imperfect Monte Carlo simulation of the momentum 
distributions for the tau decay products can also bias our 
results. The momentum distributions for charged and 
neutral pions for both data and Monte Carlo are shown 
in figure oc and figure pd. The agreement is quite good. 
Varying the slope of these distributions by ±1ct permits 
an estimate of the systematic error on A to be made. The 
effect is negligible (~ 10 -5 ). 

A possible asymmetry in the £ distribution due to 
background events has been estimated with Monte Carlo 
and found to be negligible. 
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Figure 5: The distribution for the CP-sensitive observable £ for 
each of the three choices for the scalar form factor: a) f s = 1, b) 
f s = BW[a (980)], and c) f s = BW[a (1450)]. The squares are the 
data and the line is SM Monte Carlo. 



5.2 Systematics 

Care is taken to verify that the detector does not create 
an artificial £ asymmetry. Systematic effects due to a 



6 Search for Scalar Mediated r-Decays 

We also search in a model independent way and without 
using the optimal observable £ for the presence of scalar 
mediated tau decays with the same final state as before. 




Figure 6: The ■k momentum distributions for data are shown in 
a) while the ratio is shown in b). The distributions for data and 
Monte Carlo are shown for charged pions c) and for neutral pions 
d). 
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Figure 7: The pseudo-helicity distribution for charged taus is shown 
in a) while their difference is shown in b). 



The term in Eq. ra containing Im{K) changes sign 
for taus of opposite signs. Using the angle 9*, the differ- 
ence between pseudo-helicity distributions for oppositely 
charged taus is given by 



dN(r~) dN(r+) 



d cos 9* 



d cost 



2c 2 im(A) cos 9* 



(6) 



The presence of a cos 9* term in such a distribution im- 
plies CP violation. 

The data sample used for this search is the same as 
before. The pseudo-helicity distribution for charged taus 
is shown in figure Ha. The structure in the distributions 
is due to the variation in detection efficiency as a function 
of 7r momentum and ir° energy. To measure or constrain 
7m(A) we need to determine the proportionality constant 
C2 . This is done with a procedure similar to the one used 
for the 3HDM case. 

The difference in pseudo-helicity distributions is 
shown in figure 0b. The slope extracted from the data 
distribution is C2-ZVn(A) = (4.2 ±3.6) x 1CP 4 , consistent 
with zero. This leads to a constraint on Im(A): 



0.033 < Jm(A) < 0.089 @ 90% CL 



(7) 



Note that this result is consistent with, but less restric- 
tive than our conservative limit shown in the first line 
of table |[ This is due to the use of a non-optimal CP- 
sensitive variable. 



In the SM, the helicity angle 9 V7r o, defined as the angle 
between the direction of the charges pion in the 7T7r° rest 
frame and the direction of the irir system in the r rest 
frame, has a distribution appropriate to the exchange of 
a vector particle: 



dN 



oc a + b cos 



(4) 



d cos 9^° 

Including a non-SM scalar exchange leads to a 
cos 0^° term in the distribution, corresponding to S- 
P wave interference that is proportional to the non-SM 
complex coupling A: 



dN 



d cos6L 



oc a+(ciRe(A)+C2lm(A)) cos #„.„.(> +& cos 2 0„ 



(5) 

Observation of a term proportional to cos 9^° would im- 
ply CP violation. However, the precise direction of the 
taus is unknown, due to the presence of a final state neu- 
trinos, so the distribution in Eq. o cannot be formed. 
Instead, we use the 'pseudo-helicity' angle 9* , defined as 
the angle between the charged pion in the irir rest frame 
and the direction of the irir system in the lab frame. This 
doesn't change the form of Eq. [|, just the coefficients. 
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